Quasiblack holes with pressure: relativistic charged spheres as the frozen stars 
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In general relativity coupled to Maxwell's electromagnetism and charged matter, when the 
gravitational potential W 2 and the electric potential field <j> obey a relation of the form W 2 = 
a (—£</> + b) 2 + c, where a, b and c are arbitrary constants, and e = ±1 (the speed of light c and 
Newton's constant G are put to one), a class of very interesting electrically charged systems with 
pressure arises. We call the relation above between W and <j>, the Weyl-Guilfoyle relation, and it gen- 
eralizes the usual Weyl relation, for which a = 1. For both, Weyl and Weyl-Guilfoyle relations, the 
electrically charged fluid, if present, may have nonzero pressure. Fluids obeying the Weyl-Guilfoyle 
relation are called Weyl-Guilfoyle fluids. These fluids, under the assumption of spherical symmetry, 
exhibit solutions which can be matched to the electrovacuum Reissner-Nordstrom spacetime to yield 
global asymptotically flat cold charged stars. We show that a particular spherically symmetric class 
of stars found by Guilfoyle has a well-behaved limit which corresponds to an extremal Reissner- 
Nordstrom quasiblack hole with pressure, i.e., in which the fluid inside the quasihorizon has electric 
charge and pressure, and the geometry outside the quasihorizon is given by the extremal Reissner- 
Nordstrom metric. The main physical properties of such charged stars and quasiblack holes with 
pressure are analyzed. An important development provided by these stars and quasiblack holes is 
that without pressure the solutions, Majumdar-Papapetrou solutions, are unstable to kinetic per- 
turbations. Solutions with pressure may avoid this instability. If stable, these cold quasiblack holes 
with pressure, i.e., these compact relativistic charged spheres, are really frozen stars. 

PACS numbers: 04.40. Nr,04.20.Jb, 04.70.Bw 



I. INTRODUCTION 



Frozen star was a name advocated in the 1960s by 
Zel'dovich [l[ to give to an object which is now called, 
after Wheeler's suggestion in 1968, a black hole, see e.g., 
0, see also [|| for the historical evolution of the concept. 
As we will see, a frozen star is actually a quasiblack hole. 

Quasiblack holes are objects on the verge of becoming 
black holes but avoid it, their boundary approaches their 
own gravitational radius as closely as one likes. They ap- 
peared in Einstein-Maxwell systems with special matter, 
as had been found not so explicitly in [J (see also Q) 
and then thoroughly discussed in [Q-Q . They also arise, 



and indeed were first discussed as such, in the context of 
Einstein — Yang-Mills-Higgs systems fiol . HH ] . These ob- 
jects have well-defined properties [P3 - tla |; for instance, in 
general, regular quasiblack holes are extremal. Rotating 
quasiblack holes are also known, see [l7| with hindsight 
and [H, HH for the properties of such objects. 

There are several ways of solving Einstein's equations 
when the gravitational field is coupled to matter. Since 
we will focus on Einstein-Maxwell static systems with 
matter we can perhaps distinguish two ways. One we can 
call Weyl's way, the other is the Tolman-Oppenheimer- 
Volkoff (TOV) way. 

Let us deal with Weyl's way now. Weyl [20| . while 
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studying stationary electric fields in vacuum Einstein- 
Maxwell theory, perceived that it is interesting to con- 
sider a functional relation between the metric potential 
g t t = W 2 (x z ) and the electric potential 4>(x z ) (where x 1 
represent the spatial coordinates, i = 1, 2, 3) given by the 
ansatz W — W(<p). By assuming the system is vacuum 
and axisymmetric, Weyl found that such a relation must 
be quadratic in <f>. One can go beyond vacuum solutions, 
and consider fluids which obey the Weyl relation, obtain- 
ing their properties. These fluid systems were explored 
later by Majumdar and Papapetrou [U, [22j who in go- 
ing beyond vacuum found equilibrium configurations for 
charged extremal matter where the electric repulsion is 
balanced by the gravitational attraction, and moreover 
for this particular case they showed there is a perfect 
square relationship between W and (f>, see also [23| for 
the extension of these results to d dimensions. An inter- 
esting development on Weyl's work was performed much 
later by Guilfoyle [24[ who considered charged fluid dis- 
tributions with the hypothesis that the functional rela- 
tion between the gravitational and the electric poten- 
tial, W = W(<j>) is given by W 2 = a (— e<f) + b) 2 + c, 
where a, b and c are arbitrary constants, and e = ±1 
(the speed of light c and Newton's constant G are put to 
one). This relation generalizes the usual Weyl relation, 
for which a = 1, and it allows a further beautiful rela- 
tionship between the various field and matter quantities 
[25j . which in turn generalizes the Gautreau and Hoff- 
man results [26| for fluids obeying a pure Weyl relation. 
For both, Weyl and Weyl- Guilfoyle relations, the electri- 
cally charged fluid, if present, may have nonzero pressure, 
and this turns out to be important in our context. Flu- 
ids obeying the Weyl-Guilfoyle relation are called Weyl- 
Guilfoyle fluids. 

Up to now we have only mentioned properties of the 
fluid itself, be it a Weyl-Guilfoyle fluid or, in the par- 
ticular case of zero pressure, a Majumdar-Papapetrou 
fluid. But matter solutions can be matched into vac- 
uum asymptotically flat solutions, yielding star solutions 
which, besides having the local properties associated to 
the fluid itself, have global properties for the spacetime as 
a whole. Thus, Majumdar-Papapetrou matter solutions, 
when joined to vacuum solutions yield star solutions, the 
Bonnor stars 4, 5]. These Bonnor stars when sufficiently 
compact show quasiblack hole behavior . These stars 
and quasiblack holes have no matter pressure, only elec- 
tromagnetic pressure. But now, Weyl-Guilfoyle fluids 
have, besides electromagnetic pressure, matter pressure. 
Under the assumption of spherical symmetry, Guilfoyle 
24] exhibited solutions which can be matched to the elec- 
trovacuum Reissner-Nordstrom spacetime to yield global 
asymptotically flat stars, i.e., charged stars with pressure. 

Here we explore one particular class of those spheri- 
cally symmetric cold charged fluid stars and show that 
these stars display quasiblack hole behavior, i.e., the mat- 
ter boundary approaches its own gravitational radius (or 
horizon) in a well-behaved manner. Although the cold 



charged stars have a nonextremal outer metric, the qua- 
siblack hole regime is extremal always. We analyze in 
which cases the energy conditions are obeyed, and also 
study a subclass for which the speed of sound in the 
matter is less than the speed of light. So a quasiblack 
hole with pressure obeying the necessary physical re- 
quirements is presented here. Quasiblack holes purely 
supported by electrical charge have been known [6|-|8( 
(see also [J, Q with hindsight, and @ for further results 
and references). The presence of pressure in quasiblack 
hole solutions is important, since it tends to stabilize 
the system. Indeed, an important development provided 
by these solutions is that without pressure the matter, 
Majumdar-Papapetrou matter, is unstable to kinetic per- 
turbations. So the whole solution is unstable to these per- 
turbations. Charged solutions with pressure, and along 
with them quasiblack holes with pressure, avoid this in- 
stability. Thus, if indeed stable, these compact relativis- 
tic charged spheres in the form of cold quasiblack holes 
with pressure, are really frozen stars (HQ- As a black 
hole, a quasiblack hole freezes to observers outside, but 
unlike a black hole, at the horizon limit, the star is still 
intact, it does not collapse. In the quasiblack hole con- 
text, unlike in the black hole case, the star is not irrel- 
evant. Thus the name frozen star is appropriate. From 
a well-motivated chain of works, applying Weyl's leading 
method of solving Einstein's equations we thus arrive at 
the concept of frozen stars, solutions which exhibit the 
behavior of quasiblack holes with pressure. Of course, 
quasiblack holes considered as frozen stars are more akin 
to the dark stars of Michell and Laplace (see [1]) than 
black holes themselves. 

The other way to solve Einstein's equations with 
charged matter with pressure is through the integration 
of the TOV equation for the gradient of the pressure, in 
conjunction with the other equations for the other fields. 
There are many works that use this method, we only 
quote the most relevant to our work. In an important 
work, de Felice and collaborators [13, HH found relativis- 
ts charged sphere solutions with pressure by the TOV 
method. Moreover, they took the limit to the black hole 
regime, and through numerical integration still found in 
this limit a solution. With hindsight this solution is an 
extremal quasiblack hole with pressure solution. The 
equation of state for the fluid considered is one for an 
incompressible fluid, with constant energy density, and 
in such a case the speed of sound is infinite, which might 
be considered a drawback. Bonnor [5j took notice of this 
solution and made comments and comparisons showing 
that his previous solution Q had indeed the same grav- 
itational overall properties. The search of quasiblack 
hole solutions through the TOV method is nontrivial. In 
fact, c harg ed star solutions with pressure were explored 
in [29l . |30( without finding the quasiblack hole regime. 
There are yet other ways to solve Einstein's equations. 
For instance in [ToL fll| an improved method, analytical 
and numerical, was devised to solve the Einstein — Yang- 
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Mills-Higgs system of equations. The quasiblack holes 
that arose from this system are indeed also extremal qua- 
siblack holes with pressure since the Yang-Mills-Higgs 
system has an effective built-inn pressure on its own. 
The generic properties of all types of extremal quasiblack 
holes with pressure are discussed in [l6j . 

The existence or not of solutions for compact objects in 
general relativity is connected with the Buchdahl limit. 
Buchdahl [31] found for perfect fluids that if the radius 
ro and the mass m of the star is such that m > | ro , 
then there is no equilibrium, see also [32j for more on 
this limit. For instance the Schwarzschild interior solu- 
tion with constant energy density for the matter obeys 
this inequality. If one adds electric charge to the system 
the Buchdahl limit is changed, the mass to radius of the 
star ratio ™ can certainly increase. Indeed the existence 
of quasiblack holes shows that the radius of an electrically 
charged star can approach its own gravitational radius, 
m = ?*o. Buchdahl limits for charged stars have been 
worked out in [33M36| and in particular [37[ gives a sharp 

limit, m > ^-^j - + \J if + , where q is the total 

charge of the star. For q = ro, this result admits the 
extremal case m = ro which, as we shall see below, cor- 
responds to the quasiblack hole limit satisfied by the par- 
ticular solutions studied here. Incidentally, the a —> oo 
limit of the solutions we consider yield the Schwarzschild 
interior solutions. These uncharged Schwarzschild inte- 
rior solutions contain no quasiblack hole, of course. 

Another related issue is concerned with regular black 
holes. Regular black holes are black holes devoid of sin- 
gularities. There are two types of regular black holes. In 
one type there is a magnetic core with an event horizon 
which joins into a magnetically charged vacuum solution 
different from the Reissner-Nordstrdm solution [38l - f40j 
(see also [4l[ for electrically charged regular solutions in 
nonminimal theories). In the other type, the bulk in- 
side the horizon is formed of a portion of the de Sitter 
space which joins into a Schwarzschild vacuum solution 
[43 - |47| . Curiously, the a — >• oo limit of the solutions we 
consider here yields a branch of solutions, other than the 
Schwarzschild interior solutions, which are regular elec- 
trically charged black holes. 

The present paper is organized as follows. In Sec. |TT] 
the basic equations describing a charged fluid of Weyl- 
Guilfoyle type are written. A particular spherically sym- 
metric solution to the equations representing relativistic 
charged stars, found by Guilfoyle, is shown in Sec. IIIII 
Then, Sec. IIVI is devoted to the study of the main prop- 
erties of these stars. We first review the definition of a 
quasiblack hole in Sec. lIV Al Then we obtain analytically 
the quasiblack hole limit of the solution in Sec. IIVBI 
In passing we take the a — > oo limit and show that 
there are two branches of solutions, one corresponds to 
the Schwarzschild interior solutions, the other to regular 
charged black holes. In Sec. HVCl we plot, for three dis- 



tinct cases, the relevant curves for each one of the metric 
and electric potentials, and for the fluid quantities. All 
three cases contain relativistic star solutions and quasi- 
black hole (frozen star) solutions. In Sec. [V] we conclude. 



II. WEYL-GUILFOYLE CHARGED FLUIDS: 
BASIC EQUATIONS 

The cold charged fluids considered in the present work 
are described by Einstein-Maxwell equations, which can 
be written as 

Gfxv = 8tt (T^ + £7 M „) , (I) 
VyF^ = 4tt J 11 , (2) 

where Greek indices fi, v, etc., run from to 3. G pv = 
Rfiv — \9fiuR is the Einstein tensor, with R pL , being the 
Ricci tensor, g^ v the metric tensor, and R the Ricci 
scalar. We have put both the speed of light c and the 
gravitational constant G equal to unity throughout. E pl/ 
is the electromagnetic energy-momentum tensor, which 
can be written in the form 

An = F/F vp - X -g^F piy FP a , (3) 

where the Maxwell tensor is 

Ffiu = V^X - V v Ap , (4) 

V M representing the covariant derivative, and A p the 
electromagnetic gauge field. In addition, 

= Pc U p , (5) 

is the current density, p e is the electric charge density, 
and U p is the fluid velocity. T pv is the material energy- 
momentum tensor given by 

Tpv = (p ro +p)U l _ l U u + pg pv , (6) 

where p m is the fluid matter-energy density, and p is the 
fluid pressure. 

We assume the spacetime is static and the metric 

ds 2 = g^v dx^ dx v (7) 

can be written in the form 

ds 2 = -W 2 dt 2 + hij dx i dx i , i, j = 1, 2, 3 . (8) 

The gauge field and four-velocity U p are then given 
by 

A< = -K' (9) 
U^ = -W8l- (10) 
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The spatial metric tensor hij , the metric potential B and 
the electrostatic potential <f> are functions of the spatial 
coordinates x l alone. 

The particular relativistic cold charged stars we are 
going to study here belong to a special class of systems in 
which the metric potential W and the electric potential <f> 
are functionally related through a Weyl-Guilfoyle relation 
(see Guilfoyle [24j, see also Lemos and Zanchin [25j ) 



W 2 = a (-e<]> + by + c, 



(11) 



with a, b and c being arbitrary constants and e = ±1, the 
parameter a being called the Guilfoyle parameter. Then, 
the charged pressure fluid quantities p m , p and p c satisfy 
the constraint [25[ 

abp e = e [p m + 3p+ (1 - a) p cm ] W + eap e <f>, (12) 



or 



ap c (-ect) + b) = e [p m + 3p + e (1 - a) p em ] W , (13) 

with p cm standing for the electromagnetic energy density 
defined by 



Pc 



1 (V,; 



8tt W 2 

Such matter systems we call Weyl-Guilfoyle fluids. 



(14) 



4>o being an arbitrary constant which defines the zero 
of the electric potential, and that, in asymptotically 
Reissner-Nordstrom spacetimes as we consider here, can 
be set to zero. 

Now we review the general properties of the spherically 
symmetric solutions with the boundary condition given 
above. The first integral of the only nonzero component 
of Maxwell equations ([2]) furnishes 



Q(r) - r 2 



(20) 



where the prime denotes the derivative with respect to 
the radial coordinate r and an integration constant was 
set to zero. Q(r) is the total electric charge inside the 
surface of radius r. The class of solutions we are inter- 
ested in here has c = in Eq. (|11[) . and is classified as 
class la in (24j, and so we can write (f>[r) in terms of B(r) 
as, 



e<j)(r) = b — 



B(r) 



(21) 



With this result, the amount of electric charge inside a 
spherical surface of radius r, cf. Eq. (|2U)) . is now given 
by 



Q(r) = 



-er 2 B'{r) 
2^/aMr)B(r) 



(22) 



III. SPHERICAL SOLUTIONS: GENERAL 
ANALYSIS 



We repeat in this section one class of the spherically 
symmetric solutions found by Guilfoyle [24| in order to 
study its properties. The metric (J7J instead of being 
generically written as (|5]l. is more conveniently written 
in a spherically symmetric form, namely, 



ds 2 = -B(r) dt 2 + A(r) dr 2 + r 2 dtt . 



(15) 



where r is the radial coordinate, A and B are functions 
of r only, and dfl is the metric of the unit sphere S" 2 . The 
gauge field and the four-velocity are now 



A,, 



-4>(r)8 



U, = -y/B{r) 5° 



(16) 
(17) 



The cold charged pressure fluid is bounded by a spherical 
surface of radius r = ro , and in the electrovacuum region, 
for r > ro , the metric and the electric potentials are given 
by the Reissner-Nordstrom solution 



B(r) 



1 



A(r) 
q - + 4 

r 



1 



2m 
r 



(18) 
(19) 



The continuity of the metric functions on the surface 

2?7z cj 2 

r = r yields B(r ) — 1/A(r ) = 1 1 — g-. Using 



ro 



such a boundary condition and Eq. 



o 



we get 



q_ 

rn 



'"o 



2(o- l)m 



2( a -l)-I° + 
m 



-^4a(a-l) (l 



ro 
m 



1/2 



or inverting, 



— = (1-a) — 
q r 




l + (a-l)% 



(23) 



(24) 



where the fact that Q(r = r ) = q was taken into ac- 
count. The equivalent constraints (|23|) or <\2A\ hold for 
all spherically symmetric charged pressure fluid distribu- 
tions whose boundaries are given by a spherical surface 
of radius r = ro, and whose potentials are related as in 
Eq. (|21[) . The extremal relation q = m holds when a = 1, 
in which case the relation among B and cj> is always a per- 
fect square and the ro = m limit represents a quasiblack 
hole. (In four dimensions, these Majumdar-Papapetrou 
charged fluids were studied in [2ll . |22| , and their corre- 
sponding Bonnor stars and quasiblack holes in as 
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well as in [L2T4l5ll. and in d > 4 dimensions these fluids 
were studied in [23| and their corresponding Bonnor stars 
quasiblack holes in Now, for a ^ 1, we are still con- 
sidering B = B(<p) as a perfect square (since c = 0, see 
Eq. (TTTTO but, unlike a Bonnor star, the relation q = m 
does not hold in general. However, from Eq. (f23"|) . or 
Eq. ([24"lh one sees that, for a ^ 1, there is also the possi- 
bility in which the equality mro = holds, which may 
yield back, eventually, the relation q = m. As we will 
see below, this is the quasiblack hole limit. Since the 
inequality a / 1 holds in this case, it means the corre- 
sponding quasiblack hole is made of a fluid with nonzero 
pressure. 

Guilfoyle's solutions are found under the assumptions 
that the effective energy density p c s(r) = p m (r) + 
n2 (r) . 



Pem(r) = Pm(r) 



8tt r 4 



is a constant, and that the met- 



ric potential A(r) is a particularly simple function of the 
radial coordinate. Namely [24j ]. 
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■. -l 



(25) 
(26) 



where R is a constant to be determined by the junction 
conditions of the metric at the surface r = Vq. In fact, 
by joining the function A(r) in Eq. (|26[) with the g rr 
coefficient of the exterior metric (fl"8|) . the constant R in 
Eq. (|2"5|) is found in terms of the parameters r , m, and 
q, through the equation 



1 

R? 



1 



2m 



T_ 

ro 



(27) 



are defined respectively by 



F(r) = ^\/l-^-fci 



(32) 



Q(r) 



47r / p c (r) 
h 



r dr 



l - r T 2 



d<f>{r) 
dr 



(33) 



with ki being another integration constant. The inte- 
gration constants fco and fci are determined by using the 
continuity of the metric potentials A(r) and B(r) and the 
first derivative of B(r) with respect to r at the boundary 
r — ro- The result is 



kn = 



\q\a 



2/a 



m 
(1 



q_ 

ro 



1 - 



R 2 



k R 2 - 



l-2/a 



2 - a 



I - 



R 2 



-X/a 



(34) 



.(35) 



From Eqs. (l30l) and (|33|) we get both the electric charge 
density p e and the electromagnetic energy density p om , 



87rp e (r) = 



87rp em (r) 



( v a k 2 r 2 ( 3F(r) 



2- a F 2 (r 



(2 - a) 2 F 2 (r) 



1-^-1 ,(38) 



(37) 



Since c = in the solutions we are considering (see 
Eq. (|ni)), Eq. $FZ$i with Eq. (JTTJ) turns into 



(38) 



a p c = e [p m + 3p + (1 — a) p c 



With these simplifying hypotheses, Guilfoyle [24[ was 
able to write some interesting exact solutions. We shall 
show that the class la of such solutions admits an ex- 
tremal case such that q = m = tq, which represents a 
quasiblack hole. Let us then write here this la class of 
solutions, 



B(r) 
8irp m (r) 

Q(r) 
87rp(r) 



a 



F(r) 



y+i/a 

3 a 
R 2 



2a/(a-2) 



k 2 r 2 



{2-af F 2 (r) ' 
e^/a ko r 3 



2-aF{r)' 
I a 



F 2 (r 



2k a y/l-r 2 /R 2 



(28) 
(29) 
(30) 

(31) 



2 - a F(r) 

where fco is an integration constant, and F(r) and Q{r) 



which can be interpreted as an equation of state which, as 
one can check, is obeyed by the charged fluid represented 
by the solution presented above. Furthermore, the above 
solution is valid for all a > 0, the case a = oo yielding 
the uncharged (q = 0) Schwarzschild interior solution. 

Another important quantity to determine is the speed 
of sound within the fluid. We take the usual definition 
for the speed of sound c s , 



Sp 
Sp m 



(39) 



and consider variations of the pressure p and of the energy 
density p m in terms of the radial coordinate r, i.e., 5p — 
p'(r)Sr, etc., where the prime denotes derivative with 

respect to r. Hence, we may write c 2 = — . The result is 

Pm 

the speed of sound as a function of the radial coordinate, 



b ' ' \-k 2 ^J\~r 2 lR 2 
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where we have defined ki = t. 

KqR- 1 

In the following we analyze in some detail the physical 
properties of this charged solution, exploiting, in particu- 
lar, the dependence of the metric and electric potentials, 
and of the fluid quantities on the free Guilfoyle parameter 
a. 



B. Analytical study of quasiblack holes with 
pressure (or frozen stars) 

An interesting property of the solution presented in 
Sec. Hill (cf. Eqs. pp ]) -([gT ]) ) is that all the physical quan- 
tities are well-behaved even in the quasiblack hole limit. 
In what follows we first verify that the quasiblack hole 
limit can be attained, and then we check that in such a 
limit the physical quantities related to the charged fluid 
are well-defined. 



IV. RELATIVISTIC CHARGED STARS AND 
QUASIBLACK HOLES WITH PRESSURE (OR 
FROZEN STARS) 



The metric potentials and the electric potential 



A. Definition of generic relativistic charged stars 
and definition of quasiblack holes with pressure (or 
frozen stars) 



A generic relativistic cold charged star, or sphere, is 
here defined as a smooth ball in which the gravitational, 
electromagnetic and matter fields have nonsingular be- 
havior throughout the matter, which in turn is matched 
smoothly to the Reissner-Nordstrom exterior solution. 

A quasiblack hole or frozen star can also be properly 
defined. For such a task we follow Ref. [l2j and con- 
sider the metric in the form (| L5[) with solution given in 
Eqs. (|2l)|) and (|28[1 for the interior which matches contin- 
uously to the exterior metric in the form (fT8|) . Then, con- 
sider the following properties: (a) the function 1/A(r) at- 
tains a minimum at some r* ^ 0, such that 1/A(r*) = e, 
with e « 1, this minimum being achieved either from 
both sides of r* or from r > r* alone, (b) for such a 
small but nonzero e the configuration is regular every- 
where with a nonvanishing metric potential B(r), and 
(c) in the limit e — > the potential B(r) — > for all 
r < r* . These three features define a quasiblack hole and, 
in turn, entail some nontrivial consequences: (i) there 
are 3-volume regions, rather than 2-surface regions (as 
in the black hole case), for which the redshift is infinite; 
(ii) when e — > 0, a free- falling observer finds in his own 
frame infinitely large tidal forces in the whole inner re- 
gion, although the spacetime curvature invariants remain 
perfectly regular everywhere; (iii) in the limit, outer and 
inner regions become mutually impenetrable and disjoint; 
and one can also show that (iv) for external far away ob- 
servers the spacetime is virtually indistinguishable from 
that of an extremal black hole. In addition, well-behaved 
(with no infinite surface stresses) quasiblack holes must 
be extremal. The quasiblack hole is on the verge of form- 
ing an event horizon, but it never forms one, instead, a 
quasihorizon appears. For a quasiblack hole the metric 
is well-defined and everywhere regular. Nonetheless, the 
properties that arise when e = have to be examined 
with care for each model in question. 



Let us show here that the quasiblack hole really exists 
as a special limit of the relativistic charged star solution 
presented above. According to the definition of the pre- 
ceding section, the quasiblack hole should be extremal, 
so that the mass m approaches the charge from above, 
m — > q + . In such a limit, the exterior metric (jl8p tends to 
the quasi-extremal Reissner-Nordstrom metric, for which 
the two horizons r± = m ± \J m 2 — q 2 are very close to 
each other, r + ~ m ~ r_. Moreover, there must be a 
quasihorizon r* , and then the radius of the star ro must 
coincide with r* . Hence, the quasiblack hole limit of the 
star corresponds to the limit tq — > r+ from above. We 
can then assume a relation of the form q ~ (1 — -y/e) ro, 
for a small nonnegative e. This means that the quasi- 
black hole limit corresponds to the most compact charged 
star we can have, with ro/m close to unity. In the 
present case, the mass m and the charge q are related 
by Eq. (f24|) . which implies in m/q ~ 1 + (a — l)e/2a, and 
then m/ro ~ 1 — \/s. From these results we may obtain 
the quasiblack hole limit for all other quantities. For in- 
stance, the difference m/q — q/ro which appears in the 
constants ho, ki, etc., is of the order of y/e. Taking these 
considerations into the corresponding equations we find, 
for instance, R- 2 ~ (1 - 2y/e)/rl, k ~ a 2 l a e {a - 2 ^l 2a lr\, 

B(r ) ~ £, fel ~ a (a+4)/2a £ (a-2)/2a _ fl )_ 

tion F(r) goes as k\ for all r inside the star, and then 
Eq. (|28p implies that B{r) is of the order of e for all r in 
the interval < r < ro. Moreover, we find l/A(ro) ~ e, 
satisfying the properties of a quasiblack hole as defined 
in Sec. IIV Al In the quasiblack hole limit, the metric 
potentials become 



B(r) 



A(r) 




, r <r , 
a (41) 



(42) 



—F= ) r > r 0, 
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where we have used the fact that |g| ~ m = (1 — y/e) ro 
and have written |g|/?o = l/ r o- The electric potential 
</>(r) at the quasiblack hole limit is obtained from the 
above result for B(r) and from Eq. (|21l) . 



2. The fluid quantities 



The study of the physical properties satisfied by the 
energy density and by the pressure for the present solu- 
tion was in part done in Ref. [Hj]. There are, however, 
further interesting aspects that should be considered. 

We can for instance obtain the analytical expressions 
at the quasiblack hole limit for the fluid quantities like 
mass density, charge, and pressure, 




The quasiblack hole limit of the other fluid quantities, 
namely, of the charged density p e and of the electromag- 
netic density p era are easily obtained, respectively, from 
Eq. (|38]l and from the identity p cm = Q 2 (r)/87rr 4 , and 
then we do not write them here. 

Let us now analyze the speed of sound c s within these 
rclativistic charged stars in the quasiblack hole limit. For 
that we use Eq. (|4T))) and get the limiting values of the 
constants R 2 , k , and hi, as done in Sec. IIVB II Then 
we find 



cUr) 



a + (2- 



1 



r 2 /tq 



2 - a + yfay/1 - r 2 /r, 



(46) 



As expected the speed of sound of the quasiblack hole is 
zero if a = 1. At the surface, c 2 tends to — 1 + a/ (2 — a) 
which reaches unity for a = 4/3. The function c 2 (r) at 
the center (r = 0) of the quasiblack hole is such that 
c 2 (0) is bounded to —1 from above as a tends to zero. 
In fact we see that c 2 (0) is zero for a — 1, and tends 
monotonically to —1 as a goes to zero. Hence, c s is 
undefined for all a < 1. One can further impose that the 
speed of sound is smaller than the speed of light to yield 
a further interesting class of solutions. 



3. The mass to radius relation and the a — > oo limit (the 
Schwarzschild interior and the regular black hole branches) 



The mass m, the charge q, and the radius ro of the 
relativistic stars analyzed here are related by Eq. ([23]), 
or equivalently, by Eq. (pM)) . Hence, besides the param- 
eter a, out of the three quantities m, q, and ro, we are 
left with two more free parameters. One can consider 
various possibilities. One possibility is to normalize the 
quantities in terms of the mass m. This choice is in- 
teresting, since the uncharged limit (the Schwarzschild 
interior solution) is easily obtained by taking the limit 
a — > oo. Therefore, in order to observe this limit to the 
(uncharged) Schwarzschild stars, we will in general nor- 
malize the quantities in terms of the mass m of the star. 
Another possibility is to normalize the quantities in terms 
of the charge q. This choice is also interesting since from 
it one can read directly the star's mass to radius relation, 
m x 7*0 . Of course, other possibilities could be considered. 

In order to see the difference between the two nor- 
malizations just mentioned we plot Fig. [TJ The plot on 
the left-hand side of Fig. Q] shows the behavior of q/m 
as a function of ro/m. From the curves, obtained from 
Eq. (|23"]) . we read that, for a fixed mass, the electric 
charge of the star decreases with its radius, the maxi- 
mum value being q/m = 1 for ro/m = 1, the quasiblack 
hole limit. The plot on the right-hand side of Fig. Q] 
gives the mass to charge ratio m/q for a few values of 
the parameter a as a function of the normalized radius of 
the star ro/q, whose curves are obtained from Eq. (|24|) . 
It shows that, for a fixed charge, the radius of the star 
decreases as the mass decreases. This mass to radius re- 
lation behavior is analogous to the behavior of the mass 
to radius relation in main sequence stars, and contrary 
to the mass to radius relation in white dwarfs. 

Another point to be noted is that these kind of stars 
bear a relatively large charge to mass ratio. The ratio 
m/q runs from unity, the minimum value, at the quasi- 
black hole limit, to \fa for an extremely sparse star, with 
ro — » oo. In fact, we see from Eq. (IM1) that, for fixed a, 
the maximum value of m/q is found in the limit ro — > oo , 

• • • m 
and it is given by — 

q 

Now we show that in the limit of a — > oo one recovers 
the Schwarzschild interior solution q = and p m = const, 
and also picks up another branch. Indeed, taking a — > oo 
in Eq. $Z5§ we find 



nm — - = — 

a— >oo m 2 



Because of the absolute value \ro/m — 2| one concludes, 
remarkably, there are two branches, ro/m > 2 and 
r Q /m < 2. 



( 2+ 


r»_ 2 


--) 






m 




m 
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For ro/m > 1 one finds 



lim 



0, 



(48) 



which gives the uncharged Schwarzschild interior solu- 
tion. 



For ro/m < 2 one finds 



lim 



= 2 



( 2 ~-) 

V m/ 



(49) 



This solution is weird in this context. Equation (|4TJ)l is 
equivalent to 



2m 
ro 



1- — + h=0, 



(50) 



which shows that, for this branch, when a — > oo for given 
m and q there are charged solutions in which the radius of 
the star ro is the horizon radius. For 1 < ro/m < 2 these 
are nonextremal solutions which represent noncxtremal 
regular black holes with tension matter. For ro/m — 
1 these are extremal solutions which represent extremal 
regular black holes with tension matter. Note that these 
extremal solutions can be obtained from Eq. (|4f [) when 
one takes the limit a — > oo and absorbs the e term into 
the time coordinate so that B(r) is continuous at ro- This 
a — > oo branch, representing regular black holes rather 
than quasiblack holes, should be carefully handled and 
we do not explore it further here. Regular black holes 
either with a charged core [38T - l4l| or with a de Sitter core 
[42Tj47j are known, but with charge and a de Sitter core 
together, as found here, seem to have not been explored. 



a= 1.00001 





FIG. 1: Left: the charge of the star (normalized to the mass) as a function of the radius of the star (normalized to the mass) 
for a few values of a in the interval of interest, 1 < a < 4. Right: the mass of the star (normalized to the charge) as a function 
of the radius of the star (normalized to the charge) for a few values of a in the interval of interest, 1 < a < 4. 



C. Numerical study of relativistic charged stars 
and quasiblack holes (or frozen stars): Three typical 

cases 



1 . Intervals of a 

Although valid for all a > 0, there are some special 
intervals in the domain of the parameter a for which the 
solutions can be considered more physical. R is known 
that the Schwarzschild interior solution, a = 00, having 
an incompressible fluid as matter source, yields violation 
of the dominant energy condition as well as an infinite 
speed of sound c s , i.e., a speed greater than the speed 
of light, bringing into question causality issues. Thus, 
when discussing the interesting intervals of the Guilfoyle 
parameter a, which yield what might be considered phys- 
ical solutions, it is important to take into account the en- 
ergy conditions and the behavior of the speed of sound. R 
is a straightforward task verifying that, within a certain 



range of parameters, the fluid quantities satisfy the en- 
ergy conditions for the whole star, even in the quasiblack 
hole limit. Let us then find such a interval. 

First we note that the most restrictive conditions arise 
in the quasiblack hole limit, and then we study the energy 
conditions for this kind of frozen stars. In a charged static 
fluid, besides the fluid energy density and pressure, there 
are the electromagnetic energy density p erai the radial 
electric pressure — p em and the tangential electromagnetic 
pressures p cm . It is then useful to define the effective 
energy density and pressures of the charged fluid by 

Q(r) 2 

p cff (r) ee p m (r) + p cm (r) =p m (r) + -^— r , (51) 
pl s (r) ee p(r) - p em (r) = p(r) |^ , (52) 

pl s (r) ee p(r) + p em (r) = p(r) + ^ , (53) 

where the upperscript "r" and "t" in the above defini- 
tions stand for radial and tangential pressures, respec- 
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tively. Therefore, testing the energy conditions in the 
present case leads us to check inequalities such as 

(a) PeffO) > 0, 

(b) p cS (r)+pl B (r) > 0, or equivalently, p m (r)+p(r) > 
0, 

(c) Afi(r)+Prff(r) >0, 

(d) p cS (r)+pl s (r)+2p t cS (r)>0, 

(e) Peff(r) > \p r eS (r)\, 

(f) PeffO) > bcff( r )l' or equivalently, p m (r) > \p(r)\. 

These are requirements for the weak energy condition 
(WEC) [(a),(b), and (c)], null energy condition (NEC) 
[(b) and (c)], strong energy condition (SEC) [(b), (c), 
and (d)], and dominant energy condition (DEC) [(e) and 

(f)]- 

Condition (a) is promptly verified after Eq. (|25|) . Con- 
dition (b) in the form p m (r)+p(r) > implies the charged 
fluid satisfies the inequality conditions (c) and (d). In 
addition, when the charged fluid satisfies the second in- 
equality of (f) it also satisfies (c). Therefore, it is conve- 
nient to check first the conditions on the matter quanti- 
ties p m {r) and p(r) and some relations between them. 

Starting with the matter-energy density p m (r), a nu- 
merical analysis shows that in the case ro/m = 1. 00001 
(the quasiblack hole limit) one has p m (r) > for all r if a 
is in the interval < a < 3.53523. For 3.53523 < a < 4, 
p m (r) is negative in some intervals of r inside the star, 
the corresponding (negative) minimum of /9 m (r) moves 
toward the center of the star while a increases. For a = 4, 
the minimum value of p m {r) shifts to r ~ 0.017 and be- 
comes very large (negative) but finite. For a > 4, at 
the quasiblack hole limit p m (r) behaves wildly with r. 
The extremal negative value of p-m{r) becomes arbitrar- 
ily large, moves toward the surface r — r$ as a grows, 
and eventually disappears for very large a. 

Now we turn to the behavior of the matter pressure 
p(r) at the quasiblack hole limit ro/m = 1.00001. It is 
negative for all a belonging to the interval < a < 1, 
vanishing for a = 1 in which case the matter is composed 
of charged dust. The pressure p(r) is positive in the 
interval 1 < a < 4, reaching a very large finite positive 
value at a = 4. In the interval 4 < a < oo, the behavior 
of the pressure is as wild as the energy density p m (r), 
the difference being that p(r) reaches arbitrarily large 
positive values at points where the energy density gets 
arbitrarily large negative values. 

It is also useful to compare pm{f) with p(r). For small 
a the absolute value of p(r) is smaller than p m (r). In 
particular, the central pressure is smaller than the central 
energy density, the equality p(r = 0) = p m {r = 0) being 



reached at a = (l + s/l3) /9 ~ 2.35679. For a larger 
than that value, the central pressure is always larger than 
the central energy density. A more detailed analysis valid 
for any r, not just r = 0, shows that at the quasiblack 
hole limit ro/m = 1.00001, the equality p — p m is reached 
at a ~ 2.32665, above this value one has p{r) > p m (r) 
at some interval of values of r inside the quasiblack hole 
or frozen star. Another quantity of importance involving 
Pm{?) an d p{r) together is p m (r) + p{r) for which we 
obtain that it is positive for all stars and for the frozen 
star if a is in the interval < a < 4. 

Finally, let us stress that the analysis was performed 
for ro/m = 1.00001, i.e., in the quasiblack hole limit. 
On the other hand, sufficiently sparse stars, ro/m > 1, 
satisfy all of the energy conditions for all a. 

In summary, after a careful numerical analysis of the 
quasiblack hole limit, we can state the following: 

(i) For < a < 4 conditions (b) and (c) are satisfied for 
all ro/m including the quasiblack hole. Therefore, 
the charged fluid satisfies the NEC as long as the 
Guilfoyle parameter a is in the interval < a < 4. 

(ii) The WEC requires that conditions (a) , (b) and (c) 
are satisfied. Hence it is satisfied by all of the stars 
and by the quasiblack hole or frozen star if < a < 
4. 

(iii) The SEC requires conditions (b), (c) and (d) to be 
satisfied, and so the charged fluid satisfies the SEC 
for a in the interval < a < 4. 

(iv) Conditions (e) and (f) are required by the DEC. A 
numerical analysis yields p m (r) > \p(r) \ for < a < 
2.32665 and, within this interval of a, condition (e) 
is satisfied too. Hence the DEC is satisfied in the 
interval < a < 2.33, where we put 2.32665 ~ 2.33 
from here onwards. 

Now, we analyze the speed of sound. From Eq. (|4"0")) 
we see that Cg(r) is a monotonically increasing function 
of the radial coordinate r. Hence, even if at the center 
of the star, the parameters are fixed so that the speed of 
sound is smaller than the speed of light, it may happen 
that c s (r) reaches values greater than unity at some point 
r inside the star. As the numerical analysis has shown, 
the restriction that the speed of sound is at most as large 
as the speed of light imposes the strongest bounds on the 
range of values of a. Indeed, such an imposition restricts 
the allowed values of the parameter a to be in the interval 
1 < a < 4/3. This is confirmed by the curves shown 
in Fig. [21 where we plot c s (r) as a function of a > 1, 
at the center and at the surface of the star, for several 
different values of the parameter ro/m, which represents 
the compactness of the star. The curves terminate at 
a = 1 because dp(r)/dp m (r) is negative in the interval 
< a < 1, and so the speed of sound is not defined in 
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that interval. The value a = 4/3, shown as a vertical 
dashed line in Fig. [21 is a critical value in the sense that 
for a very sparse star, i.e., in the limit ro — > oo, the 
speed of sound at the center of the star reaches unity for 
a = a c — 4/3. It is also seen from Fig. [2] that, for each 
a, the speed of sound at the center of the star strongly 



depends on the compactness parameter f3. Interestingly, 
more compact stars have smaller speeds of sound at the 
center. On the other hand, as seen from the right plots 
in Fig. [2J for a in the interval 1 < a < 4/3, the speed of 
sound close to the surface of the star is practically the 
same for all ro/m. 




1.0 1.5 2.0 2.5 




1.0 1.5 2.0 



P= 1.00001 

/?= 1.2 
- /? = 1.8 

p= 2.5 

— - ,8 = 50.0 

a 

3.0 3.5 4.0 



FIG. 2: On the left, the speed of sound at the center of the star, c s (0), as a function of the parameter a > 1, for a few 
different values of ft = ro/m. For a in the interval < a < 1, the speed of sound is not defined, since dp/dp m < 0. The value 
a c = 4/3 = 1.33 is shown as a vertical dashed line. Notice that for a very sparse star, i.e., in the limit ro — ► oo, the speed 
of sound at the center of the star reaches unity for a — a c — 4/3. For a — ► oo, the Schwarzschild interior solution, c s (0), is 
infinite at the center. On the right, the speed of sound at the surface of the star, c s (ro), as a function of the parameter a > 1. 
The value a c = 4/3 is shown as a vertical dashed line. In all the cases, i.e., for all values of /3, the speed of sound close to the 
surface of the star tends to unity at a — a c . For n-foo, the Schwarzschild interior solution, the speed of sound is infinite at 
the surface (in fact the speed of sound in the Schwarzschild interior solution is infinite throughout the whole fluid region). 



For each value of the Guilfoyle parameter a there is an 
infinity of star solutions, functions of the compactness of 
the stars themselves, i.e., functions of ro/m. The param- 
eter a is related to the pressure: for a < 1 the stars are 
supported by tension; for a = 1 the stars have no pres- 
sure, they are Bonnor stars 4-9]; and for a > 1 the stars 
are supported by pressure. In detail, a can be divided 
into the following intervals: < a < 1 yields tension 
stars; a = 1 yields charged dust stars; 1 < a < 4/3 yields 
stars and quasiblack holes which obey the energy condi- 
tions and have a speed of sound less than the speed of 
light; 4/3 < a < 2.33 yields stars and quasiblack holes 
which obey the energy conditions and have a speed of 
sound greater than the speed of light; 2.33 < a < 4 
yields stars and quasiblack holes which obey all the en- 
ergy conditions but the DEC and have a speed of sound 
greater than the speed of light; and 4 < a < oo yields 
normal stars which obey all the energy conditions but 
the DEC, but yields no quasiblack holes. In this latter 
interval of a the matter behaves as in the Schwarzschild 
solution (a — > oo), in the sense that before the gravita- 
tional radius is reached gravitational collapse ensues. 

The interval < a < 1 does not interest us here 
because it gives tension stars, the value a = 1 yields 
charged dust Bonnor stars studied previously, and the 
interval 4 < a < oo does not yield quasiblack holes and 
so again is of no interest here, although very interesting 



in other contexts. The intervals of a that yield systems 
that can be pushed into quasiblack holes is 1 < a < 4. 
So, since a = 1 has been studied in previous works, 
we will study the interval l<a<4/3, as well as the 
intervals 4/3 < a < 2.33, and 2.33 < a < 4. In order 
to see the main features of these charged stars and, in 
particular, to see the quasiblack hole limit, we have 
plotted some important curves for three typical cases 
within each interval, namely, a — 1.2, a — 1.7, and 
a = 3. We first verify that the quasiblack hole limit 
can be attained, and then we check that in such a limit 
the physical quantities related to the charged fluid are 
well-defined. 



2. The interval 1 < a < 4/3. Typical case: a = 1.2 

The metric and electromagnetic fields, and the fluid 
quantities, such as mass density p m (r), charge density 
p c {r), pressure p(r), and speed of sound c s (r), are well- 
behaved functions of the radial coordinate for 1 < a < 
4/3, even in the quasiblack hole limit. In fact, we have 
numerically analyzed each one of these functions for sev- 
eral values of ro /m as functions of the normalized radial 
coordinate rjm. Within this interval we analyze the typ- 
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ical case a = 1.2. 
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FIG. 3: The potentials and fluid quantities for the case a = 1.2. The metric potentials -B(r) and 1/A(r), and the electric 
potential <j>(r) (left panel) and the fluid quantities p m (r), p(r) and p e (r) (right panel). All quantities are plotted in terms of 
the normalized radial coordinate r/m, where m is the mass of the star, which is kept fixed, for five values of /3 — ro/m in 
each graph: j3 = 1.00001 (solid line) j3 = 1.2 (space-dashed line), /3 = 1.8 (dashed line), /3 = 2.5 (dotted line), and /3 = 50.0 
(dot-dashed line). Notice that the case ft = 1.2 is not shown for the potentials, the case ft — 3.5 being shown instead, because 
the curves of the potentials for ft — 1.2 practically coincide with the curves for ft = 1.00001. The horizontal straight dotted 
line represents the asymptotic limit of the metric potentials for large r, B(r) = 1/A(r) = 1, and cf>(r) = 0, and is plotted for 
comparison. The case /3 = 1, which gives q = m = ro, is a quasiblack hole. 



Metric and the electric potentials: The metric poten- 
tials B(r) and 1/A(r), and the electric potential (f>(r) as 
a function of the normalized radial coordinate r/m, for 
the case a = 1.2, are shown in the left panel of Fig. [3] 
The metric potentials B{r) and 1/A(r) are obtained from 
Eqs. I|28p and (f2l>)) , respectively, while 0(r) is obtained us- 
ing relation (|2"TT) with c = 0, and by choosing the constant 
b so that the potential 4>{r) is a continuous function at 
the surface r = ro- The exterior metric is the Reissner- 
Nordstrom metric, and then the curves for the potential 
B(r) and 1/A(r) tend to unity for large values of the ra- 
dial coordinate r/m. We plot the curves for each one of 
the potentials for five different values of the normalized 
radius of the star j3 = r /m, namely, /3 = 1.00001 (solid 
line), f3 = 1.8 (space-dashed line), /3 = 2.5 (dashed line), 
P = 3.5 (dotted line), and /3 = 50 (dash-dotted line). The 



case /3 = 50 represents a very sparse star and the metric 
potentials are nearly constant close to unity in the hole 
region inside the star. As a consequence, the curves of 
the potentials in this case are horizontal lines, with <f>(r) 
very close to zero. On the other hand, the electric po- 
tential is close to zero inside the star. It is also seen that 
for the quasiextremal case where ro/m = 1.00001, the 
quasiblack hole features show up. Namely, B[r) — > e for 
the whole interior region, < r < ro and 1/A(r) — > e at 
r = r . 

The fluid quantities: The mass density p m (r), the elec- 
tric charge density p e (r), and the pressure p(r) as a func- 
tion of the normalized radial coordinate r/m are shown 
in the right panel of Fig. [3l for the case a = 1.2. As 
in the case of the potential functions (see the left panel 
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of Fig. [3]) , we plot the curves in terms of the radial co- 
ordinate r/m for five different values of the normalized 
radius of the star j3 = ro/m, namely, (3 = 1.00001 (solid 
line), (3 — 1.2 (space-dashed line), (3 = 1.8 (dashed line), 
(3 = 2.5 (dotted line), and (3 = 50 (dash-dotted line). 
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quantities, as in the right panel of Fig. [3J The horizontal 
dot-dashed line is for a very sparse star, with (3 = 50. 
The radial dependence of the speed of sound is not seen 
because in this case the radius of the star is too large. 
The other horizontal line, the dashed line, indicates the 
maximum value of c s (r) close to the surface of the star 
in the quasiblack hole case, which is approximately 0.7. 
It also can be verified that for a = 4/3 such a maximum 
value is c s ~ 1, for r very close to Tq, independently 
of how compressed or how sparse the star is. In fact, 
the value of c s (ro) is the maximum value of the speed 
of sound for all the stars within the parameter space 
considered in this section. In turn, the minimum value 
of c s (r) occurs at the center of the star. 



FIG. 4: The speed of sound c s for a — 1.2, as a function of 
the normalized coordinate r/m, for five values of /3 = ro/m 
(from bottom to top: (3 = 1.00001, f3 = 1.2, /3 = 1.8, /3 = 2.5, 
and f3 = 50.0). Notice that the speed of sound is smaller than 
the speed of light throughout the star. 

The case (3 = 50 represents a very sparse star and the 
fluid quantities are very small and are nearly constant in 
the hole region inside the star, with the corresponding 
curves being the lowest horizontal lines, almost coincid- 
ing with the horizontal axes line (the bottom line) of the 
plot. In fact, such curves do not appear in the graphs. 
It is also seen that for the quasiblack hole case, repre- 
sented in the figures by the case ro/m = 1.00001, the 
fluid quantities assume the largest possible values. With 
the chosen value for the parameter a (1 < a < 4/3), 
the fluid quantities are continuous decreasing functions 
of the radial coordinate inside the star. The exterior re- 
gion is given by the Reissner-Nordstrom electrovacuum 
metric. The functions Pm(r) and p e (r) are truncated at 
r = ro, resulting in two discontinuous functions at that 
point which signal the jump into vacuum. On the other 
hand, the function p(r) is a monotonically decreasing 
function of r, starting at p(0) with the highest value and 
reaching p — at r = ro- We also see that the central 
pressure p(0) increases with the compactness parameter 
f3 of the star, analogous to the behavior of the central 
pressure in white dwarfs, and contrary to the behavior 
of the central pressure in main sequence stars. One can 
also verify numerically that the energy conditions listed 
in the previous section are satisfied within the parameter 
space considered here. 

Another final physical quantity worth of numerical 
analysis is the speed of sound. The results confirm that 
the interesting solutions, for which c s is well-defined and 
the system preserves causality, are those in which the 
parameter a belongs to the interval 1 < a < 4/3. The 
curves for c s (r) as a function of r/m for a — 1.2 are 
shown in Fig. 01 The chosen values of the normalized 
radius of the star are the same as for the other fluid 



3. The interval 4/3 < a < 2.33. Typical case: a = 1.7 

In this section we show the main features of the met- 
ric and electric functions and of the fluid quantities for a 
particular case of a in the interval 4/3<a<2. In this 
interval, all the potentials are well-behaved functions of 
the radial coordinate for all the stars. The fluid quanti- 
ties are also smooth inside the star, and satisfy the energy 
conditions (see Sec. IIVB 2[) . However, the speed of sound 
may reach values larger than the speed of light. Within 
this interval we analyze the typical case a = 1.7. 

The metric and the electric potentials: The general 
form of the curves for the metric potentials for a = 1.7 
is nearly the same as for a = 1.2 (see the left panels of 
Figs. [3] and O. This is true for all finite values of a in the 
intervals we investigated. As seen from Fig. [5] the metric 
potentials are more sensible to the radius to mass rela- 
tion, ro/m, than to the parameter a. The metric poten- 
tial B(r) in the interior region decreases slightly with a, 
indicating that the gravitational field strength increases 
with a. Similar small changes are observed also in the 
metric function A(r). The main changes, even though 
small too, are in the electric potential <p(r) which, with 
our choice of positive charge, increases with a. This in- 
dicates a noticeable change in the central electric charge 
density. 

The fluid quantities: A comparison between the right 
panels of Figs. [3J and [5] indicates that even though the 
fluid quantities strongly depend on the parameter a, for 
the values of a in the interval 4/3 < a < 2, the overall 
behavior of the density and pressure is similar to the case 
1 < a < 4/3. Namely, p m (r), Pc{r), and p(r) are mono- 
tonically decreasing functions of the radial coordinate. 
However, while the central value of the mass density 
remains almost the same, the values of the pressure and 
of the charge density increase substantially from a = 1.2 
to a = 1.7. In particular, the central pressure p(r = 0) 
increases by a factor of approximately 4. Interestingly, 
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even though the central charge density increases with 
a, the total charge inside the star decreases slowly with 
a, reaching zero as a oo for ro/m > 2 as discussed 
above. Moreover, in this interval, the energy conditions 
are satisfied by all of the stars. On the other hand, the 
speed of sound may be greater than the speed of light in 
some cases. This fact is illustrated in Fig. [5J where we 

I 
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see that the curves representing c s (r) for very compact 
stars, i.e., for j3 = r /m small, reach unity at some r 
well inside the star. It is also seen that, for sufficiently 
sparse stars, i.e., for /3 = ro/m large enough compared 
to unity, the speed of sound is larger than the speed of 
light everywhere inside the star. 
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FIG. 5: The same potentials (left panel) and the fluid quantities (right panel) as in Fig. [3] but here for the case a = 1.7. As 
above, all quantities are plotted in terms of the normalized radial coordinate r/m, for five values of /3 — ro/m in each graph, 
and the same conventions for the lines are used. Again, the case ft = 1, which gives q = m = ro, is a quasiblack hole. Notice 
that the pressure is larger than in Fig. [3] but it is smaller than the energy density. 



4- The interval 2.33 < a < 4. Typical case: a = 3.0 

Here we show the main features of the metric and elec- 
tric functions and of the fluid quantities for a particular 
case of a in the interval 2.33 < a < 4. In this interval, 
all the potentials are well-behaved functions of the ra- 
dial coordinate for all the stars. The fluid also is smooth 
inside the star, and satisfies the energy conditions (see 
Sec. IIVB 21) . However, the speed of sound reaches values 
larger than the speed of light and, for the most compact 
stars, it diverges at some point inside the star. Within 
this interval we analyze the typical case a = 3.0. 



Metric and electric potentials: The metric potentials 
B(r) and 1/A(r), and the electric potential 0(r) as a 
function of the normalized radial coordinate r/m, for the 
case a = 3, are shown in the left panel of Fig. [7J The 
general features of the curves are the same as in the case 
of the left panel of Figs. [3] and [5l for which a = 1.2 
and a = 1.7, respectively. By comparing the three cases, 
the dependence of these potentials upon the Guilfoyle 
parameter a is now more clearly seen. As a grows, the 
central values of the metric potentials B(r) and \/A{r) 
diminish by a small amount, and the electric potential 
4>{r) grows substantially when compared to the case a = 
1.2. Of course, because of the normalization used, these 
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FIG. 6: The speed of sound c s , for a = 1.7, as a function of 
the normalized coordinate r/m, for five values of /3 = ro/m 
in each graph (from bottom to top: ft = 1.00001, /3 = 1.2, 
/3 = 1.8, /3 = 2.5, and /? = 50.0). The case /3 = 1 is a 
quasiblack hole. 



changes are not observed in the quasiblack hole limit. 



The fluid quantities: The fact that the fluid quantities 
strongly depend on the parameter a is especially noticed 
for values of a in the interval 2.33 < a < 4, as seen in the 
right panel of Fig. [7J An interesting particularity is that, 
for sufficiently large a and small /3, i.e., more compact 
stars, the energy density p m (?") is not a monotonically 
decreasing function of r anymore. As seen in that figure, 
the curve of p e (r) for the quasiblack hole case oscillates, 
attaining a minimum value for some r = r m inside the 
star. This happens for all the sufficiently compact stars. 
Also worthy of note is the fact that the central value 
of p m (r) does not depend on a. Meanwhile, the pressure 
p(r) starts from a relatively high value at the center of the 
star, and decreases very fast to zero at the surface. As a 
consequence, even though the pressure is always positive 
and a monotonically decreasing function of r, the speed 
of sound increases very rapidly from the center, diverges 
at r = r rn , and becomes undefined in the region r > r m . 
For large a, the central value of the electric charge density 
also becomes much larger than the central energy density. 
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FIG. 7: The same quantities as in Figs. [3] and [5] are plotted as a function of the normalized coordinate r/m, but here for the 
case a — 3.0. The same five values of /3 = ro/m, and the same conventions are used too. As above, the case f3 = 1.2 is not 
shown in the plots for the potentials, the case j3 — 3.5 being shown instead. Notice that p m (r) does not decrease monotonically 
toward the surface r = ro- Note also that the central pressure is larger than the energy density. 



Moreover, for 2.33 < a < 4, the dominant energy con- excluding the central region of very compact stars, the 
dition is not satisfied by anyone of the stars. In fact, 
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speed of sound is greater than the speed of light for all 
the stars. This fact is illustrated in Fig. [FJ where we see 
that the curves for c s (r) are larger than unity for all r 
inside the stars. 
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FIG. 8: The speed of sound c s , for a — 3.0, as a function of 
the normalized coordinate r/m, for five values of f3 = ro/m 
in each graph (from bottom to top: /3 = 1.00001, f3 — 1.2, 
/3 = 1.8, /3 = 2.5, and f3 = 50.0). The case /3 = 1 is a 
quasiblack hole. 



V. CONCLUSIONS 



the stars are supported by tension; for a = 1 the stars 
have no pressure (they are Bonnor stars), and for a > 1 
the stars are supported by pressure. The interval of the 
free parameter a can be fixed in such a way that the 
fluid satisfies the energy conditions, and other physical 
requirements for a relativistic cold star. We have then 
studied relativistic stars within the interval 1 < a < 4. 
We have found that these cold stars show a mixed behav- 
ior, in one instance they behave as main sequence stars 
in another instance as white dwarfs. Indeed, the mass 
to radius relation of these cold stars is analogous to the 
behavior of the mass to radius relation in main sequence 
stars and contrary to the mass to radius relation in white 
dwarfs, whereas the central pressure of these cold stars 
has an analogous behavior to the central pressure of white 
dwarfs, and a contrary behavior to the central pressure 
of main sequence stars. We have also shown that, in the 
interval 1 < a < 4, the most compact configuration is 
a quasiblack hole with pressure. Thus, quasiblack holes 
without pressure (a = 1) studied previously, as well as 
quasiblack hole with pressure (1 < a < 4) studied here, 
can be found within the context of general relativity. As 
the most compact configuration is a cold star it can be 
called a frozen star. 



We have analyzed the class la of solutions provided 
by Guilfoyle (24|. Such spherically symmetric relativis- 
tic charged fluid distributions are bounded by a surface 
of radius r$. The interior region is filled with a fluid 
characterized by its mass and charge densities and by a 
nonzero pressure. The spacetime in the exterior region 
is represented by the Reissner-Nordstrom metric. These 
global solutions represent relativistic stars, i.e., relativis- 
tic cold charged spheres with pressure. Besides the mass 
m (or charge q, which are related to each other) and the 
radius of the star r^, this class of solutions is character- 
ized by another free parameter, the Guilfoyle parameter 
a. This parameter is related to the pressure: for a < 1 
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